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1 . 1 BACKGROUND 

The Navstar Global Positioning System (GPS) is a highly accurate radio- 
navigation system being developed for military and civil use by the Department 
of Defense (DoD). Due to the proposed global availability of this extremely 
accurate positioning system, GPS promises to be a major national resource for 
civil aviation and other navigation users. Of particular interest to NASA is 
a derivative of GPS, called differential GPS, which has direct applicability 
to many high-priority rotorcraft operations. Differential GPS affords 
increased levels of precision which will be essential for such rotorcraft 
applications as non-precision approach, off-shore oil operations, search and 
rescue, and oil pipeline servicing. 

Preliminary studies have been conducted by NASA to investigate differen- 
tial GPS concept mechanizations and cost, and to theoretically predict naviga- 
tion performance and the impact of degradation of the GPS C/A-code for 
national security reasons. The results of these studies show that GPS perfor- 
mance, even in the differential mode, may be inadequate (compared to FAA land- 
ing and approach navigation accuracy requirements) to support precision 
approach [1]. This is particularly true of the vertical axis accuracy of dif- 
ferential GPS, which must meet the most demanding specification. 

However, several attractive alternatives exist for improving GPS vertical 
axis performance, such as receiver "tuning** to the landing environment, 
optional selection of tracked satellites, and receiver aiding with other sen- 
sors. The objective of this effort was to investigate such techniques using 
available NASA simulation facilities, and recommend a composite system which 
meets approach and landing navigation accuracy requirements. In addition, 
flight tests were conducted with post-test differential GPS processing to 
establish a performance baseline for future flight tests of these advanced 
concepts. Results of this preliminary flight test are presented in this 
report as well. 

1.2 DIFFERENTIAL GPS CONCEPT 

The Navstar GPS is a satellite-based radionavigation system that will 
provide extremely accurate position, velocity, and time on a worldwide, con- 
tinuous basis. User receivers make ranging and Doppler measurements from 
digitally encoded L-band signals transmitted from a constellation of 18 to 21 
satellites. Four signals are required from four of the possible six to eight 
satellites in view to solve for three coordinates of position, velocity, and 
an unknown user clock bias (hence the term, pseudorange). 

User receivers employ delay lock loops and phase lock loops for tracking 
the signals and extracting the coded information to complete the triangula- 
tion-type navigation computations. In addition to residual phase and fre- 
quency errors from these tracking loops, other sources of ranging errors 
include mismodeled signal propagation delays, ephemeris errors, multipath 
errors, and intentional signal degradation imposed by the DoD for national 
security reasons. However, since the satellites are at half-synchronous alti- 
tude, the major contributor to user navigation error is non-orthogonal ranging 
signals which causes geometric dilution of precision (GDOP). 
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The sequential nature of pseudorange measurements in the receiver lends 
very well to recursive Kalman filtering, which is the usual technique for com- 
puting the navigation solution. A low-cost C/A-code tracking set will typi- 
cally employ an 8-state linearized (indirect) Kalman filter. To maximize 
efficiency and satisfy real-time constraints, the filter can be implemented 
using the upper triangular diagonal factorization of the state covariance 
matrix, using the modified Cholesky algorithm. In addition, the filter will 
need an adaptive fading memory feature to control possible filter divergence 
by increasing diagonal elements of the state covariance matrix whenever 
smoothed pseudorange measurement residuals are large. Aiding sensors can be 
useful to GPS for improving tracking bandwidths (and hence receiver noise ran- 
ging errors) and for improving the geometry of the ranging solutions. Kalman 
filter state vectors will typically be modified to estimate relevant sensor 
states. 

Differential GPS is a concept that eliminates some of the common, bias 
errors experienced by conventional GPS. Differential GPS derives its poten- 
tial from the fact that the measurement errors are highly correlated between 
different users (as well as being highly correlated in time, or autocorrela- 
tion) . By employing a second GPS receiver with comparison to truth, slowly 
varying, correlated errors can be isolated and eliminated. In addition, 
depending on the relative rates, intentional degradation of the C/A-signal may 
be eliminated by differential GPS as well [2] . A popular implementation is 
illustrated in Figure 1-1. 




Figure 1-1. Differential GPS Concept 
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1.3 ROTORCRAFT APPROACH AND LANDING 


The attractive application of GPS to the rotorcraft community is for 
remote, low visibility meteorological conditions where other reliable and pre- 
cise forms of navigation are unavailable. Specific examples include Alaskan 
inter-city navigation? off-shore oil rig navigation, approach, and landing; 
and mountainous operations where precise beacons are unavailable. Also, spe- 
cial operations such as search and rescue, police work, and fire fighting can 
benefit from enhanced navigation precision. 

For landing operations in these locations and applications, even more 
precise positioning data is required. Such precision may be available by 
employing differential GPS. Thus, a differential GPS user may be totally 
independent of visibility conditions, which for many of the locations postu- 
lated constitutes a large percentage of the operations. 

Rotorcraft approach profiles vary, but are generally flown at six to ten 
degree glideslopes. Low visibility approaches usually include long straight- 
in finals, but finals as short as one mile or less are feasible. Speeds vary 
by equipment flown, from as fast as small fixed-wing aircraft to much slower 
speeds. The controllability of speed as well as flight path is a factor that 
may have potential in relaxing landing navigation standards, currently esta- 
blished for fixed and rotor wing aircraft in general. 

However, for the present, FAA Navigation System Accuracy Standards are as 
shown in Table 1-1. As stated, the elevation accuracy requirement is more 
demanding than the lateral error. Even at Category I decision heights, the 
vertical axis requirement presents a formidable challenge to differential 
GPS. At these accuracy levels, every known error source in differential GPS 
must be exploited in full to have any hope of reducing cumulative error bud- 
gets below required levels. Furthermore, because GPS solutions often lack a 
satellite directly at zenith, the vertical eucis suffers the greatest accumula- 
tion of multiple sources of errors. 


Table 1-1. Approach and Landing Navigation Accuracy Requirements 
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II. PROBLEM DKFINITIOtl AND SCOPB 


2.1 OBJECTIVES 

The objectives of this effort were to improve the vertical axis perfor- 
mance of GPS to potentially support precision helicopter landing approach 
operations • 

Analysis of differential GPS performance have shown that while the system 
may provide adequate (with respect to current PAA standards) accuracy in the 
lateral aixis to support precision approach, further improvements are necessary 
in the vertical axis. Fortunately, investigation of GPS error sources indi- 
cates that several receiver and navigation filter design changes are possible 
which may improve vertical axis accuracy for this particular application. 

Such error models and filter components are inherent in the NASA DIFFGPS simu- 
lation, thus providing a readily available facility for verifying this hypo- 
thesis • The purpose of this effort is to determine improvements that can be 
realized in rotorcraft landing approach differential GPS vertical axis perfor- 
mance using the DIFFGPS simulation. Specific areas of investigation included 
receiver measurement processing improvements, state modeling improvements, 
adaptive covariance modeling, mission-tailored satellite selection, and inte- 
gration of external aiding sensors. 

2.2 SYSTEM SIMULATION FACILITY 

The major GPS evaluation tool used for this study was DIFFGPS. DIFFGPS 
is an analytical simulation of the Navstar GPS and its environment. In addi- 
tion to conventional GPS, the simulation models the differential implementa- 
tion of GPS, where a second, static receiver at a surveyed location is used to 
compute measxxrement corrections and augment the dyncuaic user ' s solution via 
telemetry data link. Figure 2-1 illustrates the DIFFGPS program. 

DIFFGPS was designed as a Monte Carlo simulation to provide realism of 
the environment and its non-linear features, and to provide direct relevance 
to field test programs [3] . Satellite motion is modeled by classic Keplerian 
equations. Constellation alternatives are selectable for up to 24 satel- 
lites. User motion in the User Route Planning Prograua is produced by the 
operator establishing a "route plan" by specifying either latitude/longitude/ 
altitude or range/bearing/altitude from an initial location. The simulation 
generates accelerating turns, climbs, dives and linear speed changes, automa- 
tically checking for adequate acceleration distances in the route plan. 

Particular emphasis was placed on the development of the error models, 
since they explicitly determine navigation performance. Care was taken to 
assure that empirically-backed error levels were incorporated. Dyneunics of 
the errors due to user-satellite and user-corrector geometry were faithfully 
reproduced to support representative differential GPS performance conclu- 
sions. This caused consideration and modeling of somewhat undefined relation- 
ships, such as spatial and temporal variations of ionospheric and tropospheric 
delays. Magnitudes of the errors may be adjusted so that the relative compar- 
ison of different filters should produce very robust conclusions based on the 
error model fidelity used. This last feature, reliability of the comparative 
performance of the filters, was the major design methodology driver in devel- 
oping the various error models. 
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GPS navigation algorithm analysis is the primary emphasis of the simula- 
tion. To enable an efficient and fair comparison of implementation techniques 
under the Monte Carlo driver environment ^ the simulation can run three correc- 
tor and six user navigation algorithms simultaneously during any run^ although 
only one corrector is permitted to provide the differential reference solution 
for a particular run. Therefore, five user receiver navigation algorithms can 
process identical motion, satellite geometry, and error model inputs during a 
simulation run. Due to the modular nature of the simulation, it is possible 
to substitute different sets of variations of filters into the six user and 
three corrector slots prior to the run. Of course, it is equally easy to sub- 
stitute in variations of error models or even satellite/user motion models if 
the need arises. Interface specifications are well defined for such flexibil- 
ity. 


The post-run analysis module provides the analyst extensive observation 
of relative motion, error model outputs, and Kalman filter operation. This is 
accomplished by providing menu-selectable plots, with selectable data periods 
and scales, for all parameters of interest. The operator can essentially 
retrieve all significant input, intermediate, and output variables from a 
large data file produced during the run. A high-performance plotting package 
and high-resolution graphics provide the display medium, including multiple 
plots per page and multiple curves per plot for comparative analysis. 



Figure 2-1. DIFFGPS Monte Carlo Simulation 
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III. RBCBIVKR NAVIGATIOH PILTKR KNHAWCBtBtfT 


This section describes enhancements to the receiver measurement process- 
ing and Kalman filter development. The measurement processing change involves 
continuous deltaranging and Is discussed In Section 3.1. 

The two significant improvements to the earlier modeled Z-Set type filter 
were acceleration state modeling and Improvements to adaptive covariance 
modeling. Such models Improve the fidelity of the plant model to the known 
dynamics of the vehicle, yet allow graceful boost of process noise when the 
model Is Inadequate. It Is generally true that when the plant model Is more 
precise, the process noise modeling must be responsive to preserve optimal 
filter operations In all conditions, particularly In those conditions for 
which the modeling was not designed and may be In error. Acceleration state 
modeling Improvements are discussed In Section 3.2 and Adaptive Covariance 
modeling Is discussed In Section 3.3. 

3.1 MEASUREMENT PROCESSING IMPROVEMENTS 

Since the baseline GPS receiver technology for this effort was the Magna- 
vox Z-Set, an analysis was made of the receiver measurement process to deter- 
mine possible Improvements. 

The major area for Improvement and a technique that Is being Implemented 
In modern receivers Is continuously Integrated Doppler (carrier processing). 

In the Z-Set, the effective Doppler Integration Interval Is only 1/4 of the 
pseudorange measurement Interval. This was necessary because of the delta- 
pseudorange technique of deriving Integrated Doppler. That Is, during the 1/4 
Interval Doppler Integration process, the code loop Is disabled so that the 
end-of-deltarange-lnterval pseudorange measurement Is precise to relative 
Doppler accuracies, but this technique causes the deltarange to be non- 
contiguous In time. Therefore, use of deltarange for velocity determination 
requires extrapolation which will be In error under accelerating conditions. 

In particular, the "smoothing" advantage of continuous, precise deltarange 
measuments Is lost to a great extent In this method. 

In current, multi-channel GPS receivers with sufficient tracking loop 
processing power, measurement processing can achieve full-interval deltarange 
averaging or nearly so. This deltarange Is a very precise (fraction of a car- 
rier wavelength, 19 cm, typically 1-2 cm) measure of change in pseudorange and 
therefore, has a tremendous effect on smoothing the less accurate (3-6 m, 

1 a typically) pseudorange measurements If applied correctly. Although peri- 
odic loss of carrier lock would force "reinitialization" of this smoothing 
process, the process quickly converges, so If positive coherent tracking Indi- 
cation Is present, the advantage of having contiguous deltaranging Is still 
apparent. 

To account for this tracking Implementation improvement in DIFFGPS, the 
receiver tracking model was modified to emulate full-interval deltarange 
tracking. Figure 3-1 compares the performance of an 8-state (position- 
velocity) filter tracking aircraft motion during a constant acceleration turn 
followed by a constant linear acceleration. 
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The full internal deltaranging had a dramatic effect on performance, in 
fact/ the effects of the acceleration are virtually ludetectable. This is due 
to the fact that the process noise is (emd is modeled as) nearly zero, since 
the deltarange is a very accurate measurement of the true change in velocity 
over the interval, not unlike the use of delta velocity inputs from an INS. 

The error excursions of 20-40 meters in the Z-set modeled case are typical of 
the case of partial interval doppler measurement. 


Full Interval Deltaranging 
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Figure 3-1. Comparison of Performance for 1/4-Interval 
and Full-Interval Deltarange Integration 

For acronym definition refer to glossary on page vi. 
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3.2 ACCELERATION STATE MODELING 


In general, the inclusion of acceleration states in the Kalman Filter, 
allowing for the estimation of the vehicle acceleration should provide for a 
better dynamic tracking of the vehicle motion. However, since the only direct 
observations in the measurement vector are range and range rate, any accelera- 
tion is observable only as the derivative of the rate and therefore, lags the 
velocity estimate. For this reason, acceleration state modeling will work 
well only when the vehicle acceleration is reasonably constant, allowing the 
filter to build up that state and stay there. Conversely, changes in acceler- 
ations will be a transient to this process and will cause the filter to lag in 
its response unless adequate process noise is introduced to deweight the pro- 
cess (constant acceleration) modeling. 

Of course, it is important to remember that the state modeling is related 
to the measurement and computation intervals. For example, if the measurement 
and computation intervals are 1 second, then it is important to closely model 
the dyneunics over that 1 second interval rather than be concerned with longer 
term effects. The significance of this is that, as tracking intervals get 
shorter for the same vehicle dynamics, the structure and perhaps even order of 
the state modeling becomes less important. 

The models considered in this study were the baseline posit ion- velocity 
state vector (with the contiguous deltaranging discussed in Section 3.1), con- 
stant acceleration state modeling, and a concept called constant turn-rate 
dynamics. Comparisons for a constant contripetal acceleration, level turn 
were performed using the OIFFGPS Simulation and are shown in this section with 
the various modeling discussions. The most attractive models were later 
tested with the other design improvements in the landing simulation which is 
reported in Section VI. 


3.2.1 No Acceleration Modeling; PV Filter 


This model, which emulates the baseline Z-Set Kalman filter (except for 
the addition of contiguous delta-ranging) investigated in an earlier study 
[3] , uses an eight-state filter containing the three position states and three 
velocity states, the user clock bias, euid the user clock drift rate, as 
follows: 
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where the w terms are white process noise with normal distribution, zero mean 
and power spectral density N, i.e., w^ ). This formulation models, 

or at least accounts for, acceleration^related effects in the process noise 
matrix, Q. The impact of this is to boost Q, and hence, the covariance, P, 
more than would be done for the acceleration state case which would incorpor- 
ate only jerk-related effects in the Q-matrix. With larger Q then, the posi- 
tion-velocity (PV) filter will deweight the state modeling (state extrapola- 
tion) and weight higher the measurements, resulting in greater response to 
measurement errors and noise. 

3.2.2 Constant Acceleration; PVAC Filter 


In this model, a constant acceleration is assumed in each of three 
axes. The filter consists of 11 states, including the three positions, three 
velocities, three accelerations, the clock bias and drift rate. 

The principle behind this model is that# as long as actual accelerations 
are constant, the acceleration state estimates will attain constant values. 
Constant acceleration is, in fact, a good assumption for most flight sit\ia- 
tions since power changes normally result in reasonably constant linear accel- 
erations while coordinated turns, in zero wind conditions, result in constant 
centripetal acceleration. Of course, the ECEF components of total accelera- 
tion are not constant as defined in the model, which is a disadvantage of this 
technique if modeled in the ECEF frame. 

The acceleration is modeled as a random walk, i.e., a white noise forcing 
function is used as process noise. 

The equations are: 



b » f + Wjj 

f =r -f/T + w^ 


( 2 ) 
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3.2.3 Turn Rate Dynamics; PVAT Filter 


In this model, a constant turn rate is assumed, i.e., 

a(t) » - m v(t) + w(t) 


where u is the turn-* rate, or 


( 3 ) 


(0 - v(t) X a(t) 


(4) 


v(t) 

Note that u is not a constant, but rather a variable non-linear function of 
time. 




Clearly, the use of this model will add a certain ^uaount of computational 
burden and the discussion below trades off this burden against the improved 
accuracy over simpler models. 

This model takes advantage of the assumption of coordinated turns where 
the acceleration vector is perpendicular to the instantaneous velocity vec- 
tor. The model identifies such a situation £uid adjusts the acceleration state 
accordingly. 

3.2.4 Discrete Formulation 

The continuous formulation provided in the previous sections can be sum- 
marized as: 


X = FX -t- W (5) 

The equivalent discrete formulation is obtained by 

X(k + 1) - (j> (k)X(k) + W(k) (6) 


with 


^ (k) - L"”*(sl - P)"”' 


(7) 


where is the Inverse Laplace transform and 


W(k) 


/t 


'k * ’ 




(8) 


These calculations were performed for the three models described previously 
and are included in Appendix A, along with the discrete noise covariance 
matrix. 


QH«) 


/ T) N(t) ^ T) dT 

k 


( 9 ) 


Note that all matrices are constant, except for the turning dynamics model, 
where they are functions of time (through the turning rate oj). 


I 
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3. 2* 5 Simulations 


Simulations were performed using the DIFFGPS prograun to assess perfor- 
mance for the various models during both conventional and differential GPS 
operation and for several levels of dynamics. 

Although the object of the study is vertical accuracy enhancement, the 
performance of the various filters was initially assessed using a level Con- 
stant Acceleration Turn type pattern. It was felt that accelerations were 
experienced in all axes, so that improvements laterally would imply similar 
improvements in the vertical axis. Performance in the vertical axis is des- 
cribed in later sections where an actual descending glidepath profile is simu- 
lated. The path is plotted in Figure 3-2. 

The Constant Acceleration Turn profile selected for the simulations con- 
sisted of: 

a. a straight, level, constant velocity segment at 1000 ft altitude, 100 

knots, and -80 degrees bearing (duration: 29 seconds); 

b. a level turn, constant velocity segment at 1000 ft altitude, 100 

knots to a bearing of +10 degrees (duration: 11 seconds); 

c. a straight, level, constant velocity segment at 1000 ft altitude, 100 

knots, and +10 degrees bearing (duration: 22 seconds); 

d. a straight, level, accelerating segment at 1000 ft altitude, +10 

degrees bearing, and final velocity of 155 knots (duration: 9 

seconds ) ; 

e. a straight, level, constant velocity segment at 1000 ft altitude, 155 

knots, and +10 degrees bearing (duration: 23 seconds). 

3.2.5. 1 PV Filter 

The PV filter performance is presented in Figures 3-3 and 3-4 for this 
type of process noise matrices. 

The baseline Z-Set model used a diagonal process noise matrix. Ignoring 
cross terms. While this formulation is technically incorrect, it will not 
necessarily result in poorer performance due to the reduced order state vec- 
tor, non-linearities, and measurement characteristics. In particular, Eller 
[4] pointed out in his study that the high precision of the deltarange meas- 
urements compared with pseudorange measurements (factor of 100 or 200-to-1) 
encourages decoupling of the position states form velocity (and acceleration, 
if modeled) states in the process noise matrix. 

In Figure 3-3, a diagonal process noise matrix is used, while in Figure 
3-4, the full process noise matrix is used. 
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Figure 2 - 2 . Constant Acceleration Turn Flight Profile 
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Figure 3-3a. PV Filter Performance, X-Axis Diagonal Q-Matrix 
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Figure 3-3b. PV Filter Performance, Y-Axis Diagonal Q-Matrix 
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Figure 3-3c. PV Filter Performance, Z-Axis Diagonal Q-Matrix 
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Figure 3-4a. PV Filter Performance, X-Axis Full Q-Matrix 
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Figure 3-4b. PV Filter Performance, Y-Axis Full Q-Matrix 
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Figure 3-4c. PV Filter Performance, Z-Axis Full Q-Matrix 


43-3020-18 



A comparison between the two cases indicates little difference over the entire 
path in both the conventional and differential modes. RMS and Standard devia- 
tions are listed in Table 3-2 for all three axes. 

Table 3-1. PV Filter Performance (Meters)/ 

Constant Acceleration Turn Pattern 


VALUES IN METERS 


FULL Q 
RMS 

MATRIX 

STANDARD 

DEVIATION 

DIAGONAL Q 
RMS 

MATRIX 

STANDARD 

DEVIATION 

CONVENTIONAL GPS 

x 

13.8 

4.0 

13.8 

5.7 


y 

30.6 

8.5 

22.5 

8.8 


z 

7.3 

4.3 

7.4 

4.3 

DIFFERENTIAL GPS 

X 

4.5 

4.5 

6.0 

. 

in 


y 

9.0 

7.1 

11.3 

2.0 


z 

3.4 

2.8 

3.3 

2.8 


The table shows that the differences between the two filters are not 
large and that they are not favoring one over the other. This is true for 
both the conventional and differential modes. During the filter timing/ it 
becaune clear, however, that one had to be very careful in the case of the full 
Q matrix since the presence of the off-diagonal terms could lead to an ill- 
defined Q matrix (i.e., non-positive definite). The data presented in Table 
3-1 is for the full Constant Acceleration Turn pattern. Table 3-2 presents 
the data if only the turn segment is considered and, therefore, shows better 
how each filter tracks during the higher dynamics. Here, too, differences are 
not clearly in favor of one model or the other, although the standard devia- 
tions are generally smaller for the full Q matrix filter. Overall, though, 
conclusions are difficult to draw with the PV filter as to the advisability to 
use a full or diagonal Q matrix. 
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Table 3-2. PV Filter Performance, Turning Segment of 
Constant Acceleration Turn Pattern 

FULL Q MATRIX DIAGONAL Q MATRIX 


VALUES IN METERS 


RMS 

STANDARD 

DEVIATION 

RMS 

STANDARD 

DEVIATION 

CONVENTIONAL 

GPS 

X 

10.1 

4.6 

9.3 

7.2 



y 

26.8 

6.8 

30.6 

10.1 



z 

2.1 

1.1 

2.4 

1.3 

DIFFERENTIAL 

GPS 

X 

5.9 

4.8 

9.6 

7.3 



y 

14.3 

6.9 

18.5 

10.2 



z 

5.1 

0.7 

4.2 

2.5 


3. 2. 5. 2 PVAC Filter 

Figure 3-5 presents the PVAC filter performance for the Constant Acceler- 
ation Turn pattern and using a full process noise matrix. In Figure 3-6, the 
performance is presented for the diagonal process noise matrix where the posi- 
tion states are decoupled from the velocity emd acceleration states, leaving 
cross terms between velocity and acceleration only. Table 3-3 shows the RMS 
and standard deviations for all axes and for both the conventional and differ- 
ential modes. As for the PV filter, the performance for both filters are not 
significantly different. Table 3-4 shows the performance during the turning 
segment only, and, here too, the differences are minor and do not advantage 
one implementation over the other. 

Table 3-3. PVAC Filter Performance (Meters), 

Constant Acceleration Turn Pattern 

FULL Q MATRIX DIAGONAL Q MATRIX 





RMS 

STANDARD 

DEVIATION 

RMS 

STANDARD 

DEVIATION 

CONVENTIONAL 

GPS 

X 

14.3 

2.8 

14.6 

2.5 



y 

17.8 

4.8 

17.6 

5.2 



z 

7.2 

4.4 

7.5 

4.5 

DIFFERENTIAL 

GPS 

X 

3.8 

3.7 

3.8 

3.7 



y 

6.3 

5.3 

6.6 

6.1 



z 

3.5 

2.9 

3.7 

3.0 
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Table 3-4. PVAC Filter Performance (Meters), 

Turning Segment of Constant Acceleration Turn Pattern 

FUU. Q MATRIX DIAGONAL Q MATRIX 




RMS 

STANDARD 

DEVIATION 

RMS 

STANDARD 

DEVIATION 

CONVENTIONAL GPS 

X 

15.2 

2.3 

9.2 

4.6 


y 

22.7 

2.7 

26.4 

2.8 


z 

1.7 

1.4 

3.0 

0.15 

DIFFERENTIAL GPS 

X 

4.3 

2.8 

6.1 

4.8 


y 

9.8 

2.7 

14.1 

2.8 


z 

5.5 

1.3 

4.9 

0.16 
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Figure 3-5a. PVAC Filter Performance^ X-Axis Full Q-Matrix 
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Figure 3-5b. PVAC Filter Performance, Y-Axis Full Q-Matrix 
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Figure 3-5c. PVAC Filter Performance, Z-Axis Full Q-Matrix 
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Figure 3-€a. PVAC Filter Performance, X-Axis Diagonal Q-Matrix 
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Figure 3- 6b. PVAC Filter Performance, Y-Axis Diagonal Q-Matrix 
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Figure 3-6c. PVAC Filter Performance, Z-Axis Diagonal Q-Matrlx 
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Comparing the PVAC filter performance with the PV filter performance, 
based on the full Constant Acceleration Turn pattern, performance is usually 
slightly better for the PVAC filter, although the impact of any major improve- 
ment achieved during the turning segment would be somewhat reduced due to the 
averaging with the straight segments* Comparisons of the turning segments 
alone (Tables 3-2 and 3-4) show a clear advantage (factor two or three) in 
favor of the PVAC filter for the standard deviation. This is not as clear for 
the RMS values, however. This can be explained by realizing that during peri- 
ods of poor tracking by the filter, the errors can actually grow in a way to 
reduce the absolute value of the mean over the path and therefore reduce the 
RMS, providing therefore a false indication of improved performance. The 
standard deviation, however, indicates the fluctuations in the error and is a 
better indication of the filter performance. 

3. 2. 5.3 PVAT Filter 

As was mentioned in Section 3.2.3, the turn rate is a function of time. 
Since the process noise matrix is a function of the turn rate, it, too, is 
therefore a function of time. Performance was evaluated using both a time 
variad^le and a constant process noise matrix* The diagonal noise matrix was 
used, i.e*, the position states are decoupled from both the velocity and 
acceleration states, learning cross terms between velocity auid acceleration 
only. Figure 3-7 shows the performance for the Constant Acceleration Turn 
pattern with the variable process noise matrix. Performance with the constant 
process noise matrix is presented in Figure 3-8. 

In Table 3-5, the RMS and standard deviation of the errors are provided 
for all axes and for both the conventional and differential modes, clearly, 
the effects of the variable Q matrix are insignificant in the filter perfor- 
mance. Table 3-6 shows the performance with both the constant and variable Q 
matrices during the turning segments of the Constant Acceleration Turn pat- 
tern. Here, too, the effects are very small. Time variable process noise 
matrices are therefore not necessary and certainly not worth the additional 
processing load* 


Table 3-5. PVAT Filter Performance (Meters), 
Constant Acceleration Turn Pattern 


Diagonal Q Matrix 

VARIABLE Q MATRIX 

STANDARD 
RMS DEVIATION 


CONSTANT Q MATRIX 

STANDARD 
RMS DEVIATION 


CONVENTIONAL 

GPS 

X 

14.5 

2.4 

14.5 

2.5 



y 

17.5 

5.0 

17.4 

5.2 



z 

7.5 

4.5 

7.5 

4.5 

DIFFERENTIAL 

GPS 

X 

3.8 

3.7 

3.8 

3.7 



y 

6.6 

6.1 

6.5 

6.0 



z 

3.7 

2.4 

3.6 

2.3 
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Figure 3-7a. PVAT Filter Performance, X>Axis Ulagonai y-natrix 
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Figure 3-7b. PVAT Filter Performance, Y-Axis Diagonal Q-Matrix 


43 - 3020-30 



Position Error (m 

la.D -IO.D -3.0 0.0 ! 


LTP CocrdLrusL* Sustain DDia 
18 Sat. cpheoerts Date s I JfiN 83 



TIME ISEC) 


NflV 2(+) PVflT C i -0.S99BE+0J, 0,4483E+01, 0.7487E->0] 

NflV Z(-^) PVflf 0 » 0.2H4E«-0I, 0.2S93E+01, 0.3680E^-0] 


For acronym definition refer to glossary on page vi. 


Figure 3-7c. PVAT Filter Performance, Z-Axis Diagonal Q-Matrix 
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Figure 3-8a. PVAT Filter Performance, X-Axis Constant Diagonal Q-Matrix 
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Figure 3-8b. PVAT Filter Performance, Y-Axis Constant Diagonal Q-Matrix 
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Figure 3-8c. PVAT Filter Performance, Z-Axis Constant Diagonal Q-Matrix 
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Table 3-6. PVAT Filter Performance (Meters) , 
Turning Segment of Constant Acceleration Turn Pattern 


Diagonal Q Matrix 

VARIABLE Q MATRIX CONSTANT Q MATRIX 




RMS 

STANDARD 

DEVIATION 

RMS 

STANDARD 

DEVIATION 

CONVENTIONAL GPS 

X 

13.5 

1.4 

13.6 

1.7 


y 

20.6 

4. 1 

20.2 

4.7 


z 

2.0 

1.2 

2.1 

1.3 

DIFFERENTIAL GPS 

X 

1.6 

1.5 

1.9 

1.8 


y 

6.4 

4.8 

7.8 

4.7 


z 

5.4 

0.7 

5.2 

0.9 

Con^arlng the 

PVAT 

filter performance with the 

PV and PVAC 

filters, very 

little difference 

can be 

noted for 

the full Constant 

Acceleration Turn pat- 


tern. For the Turning segment only, PVAT filter standard deviation is much 
smaller (by a factor of 2 or 3) than that obtained with the PV filter. Com** 
parison with the PVAC filter does not show, however, a clear cut advantage for 
either the PVAT or PVAC filter. At the RMS level, however, the PVAT filter 
seems to perform better although this could be incidental to this particular 
run only. 


3.2 .5.4 Conclusions for Constant Acceleration Turn Pattern 

Comparison of the PV, PVAC, and PVAT showed a clear improvement obtained 
by modeling the acceleration. The performance of the PVAC filter, which 
models the acceleration on a constant use, however, not significantly differ* 
ent than that of the PVAT filter which models the acceleration as centripe- 
tal. The performance obtained by using the full theoretical process noise 
matrix was not significantly different than that obtained using an diagonal Q 
matrix, in which position states are decoupled from both velocity and acceler- 
ation states. 

The PVAT filter performance using a time varying process noise matrix did 
not improve much over that obtained with a constant process noise matrix. It 
appears, therefore, that the PVAC filter seems best suited overall for a Con- 
stant Acceleration Turn type pattern, considering both performance and com- 
plexity factors. 

The PV, PVAC, and PVAT filters were further tested using a landing type 
pattern. Results are discussed in Section 6. 

3.3 ADAPTIVE COVARIANCE FEATURES 

In the PVAT filter, continuous centripetal acceleration is modeled. How- 
ever, since dynamics is rarely a continuous centripetal acceleration, this 
modeling could, in some cases, be worse than for a simpler filter, assuming 
either constant (PVAT) or just noisy (PV) acceleration. In that case, an 
increase in the process noise would be required to offset the modeling 
errors. High dynamics transitions (i.e., large changes in turning 
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rate, o>) should be accompanied by increased in the process noise, Q. This was 
simulated here in the following way: 

At t - 0, Q - 

At t « t + 1, 

If [(AdJ)^ greater than 10 
then t ■ 0 

Q - Q [1 + a(Ao))2 ^ 

o 

If [t«greater than. 5] 
then Q - Qj, 

Basically, if^CAcu) exceeds a threshold, the process noise is multiplied 
by [1 + a(A(u) ] and then allowed to reduce back to in five seconds, unless 
a new larger Aw occurs first, in which case, the clock gets reset to zero, Q 
gets reset to a higher value and reduces back to in five seconds* The term 
0( is a multiplicative factor that can be set by the simulation operator* 

The Constant Acceleration Turn pattern, used in the previous tests, is 
used here* The PVAT filter was run with a constant process noise (diagonal) 
matrix and then with an adaptive process noise (diagonal) matrix* The multi- 
plicative factor "a" was set to 50,000* 

In Section 3 *2* 5 *3, the performance of the PVAT filter was studied and 
plots were provided* Little difference was noted when the adaptive scheme was 
tried using the same process noise matrix, even with very large multiplicative 
factors* This would tend to indicate that the filter does not rely much on 
the model so that increased process noise is of little or no effect* 

To determine if adaptive schemes could provide improved performance, the 
PVAT filter was run with a reduced process noise matrix* The performance is 
presented in Figure 3-9. Clearly, the filter experiences difficulty during 
the turn segment* Figure 3-10 presents the performance of the PVAT filter 
using the adaptive process noise matrix with the same initial value as that 
used for the constant case, indicating a 50% reduction in the error due to the 
turn * 


Table 3-7 shows the RMS and standard deviations for all axes and for both 
the conventional and differential modes for each of the PVAT filters discussed 
above* A slight improvement is noted in one direction while some degradation 
is observed in the vertical direction when using the adaptive scheme* The 
figure shows, however, that the adaptive filter tracks better in the turn, 
even if the RMS error turns out to be larger ( due to biases ) . 
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Figure 3-9a. PVAT Filter Performance, X*Axis Reduced, Constant, 

Diagonal Q-Matrix 
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Figure 3-9b. PVAT Filter Performance, Y-Axis Reduced, Constant, 

Diagonal Q-Matrix 


43-3020-38 



Position Error (m.) 

ia.o-12.3 -10.0 -y.s -aj - 2.9 0.0 2.5 9.0 r.s in.o 


LTP CoordLnaL* Sustmn OcrLa 
18 Sot. Epheaerls Oota i 1 JflN 83 


CU-ROUfiHR COUSVY 
a - NflV Z(+) PVfll C 

♦ -Nflvz(-») pvfli a 



I'a 

1 — 

ao 

so 

TIME 

-P 1 1 

00 70 80 

(SCO 

ob 

lOQ 

NflV 

NAV 


pvflr c 
pvnr 0 

t -0.5274E+0J, 
» 0.246BE*01, 

0.3CQ2E<^01, 0.6078E-Kn 
0.l933B-*-0I, 0.3135E-^03 




For acronym definition refer to glossary on page vi. 


Figure 3-9c. PVAT Filter Performance, Z-Axis Reduced, Constant, 

Diagonal Q-Matrix 
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PVAT Filter Performance, X-Axis Adaptive 
Process Noise Q-Matrix 
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Figure 3-1 Ob. PVAT Filter Performance, Y-Axis Adaptive 
Process Noise Q-Matrix 
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Figure 3-1 Oc. PVAT Filter Performance, Z-Axis Adaptive 
Process Noise Q-Matrix 
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Table 3-7. PVAT Filter Performance (Meters) 


Adaptive Mode 


Diagonal Q Matrix 





CONSTANT Q MATRIX 

ADAPTIVE Q MATRIX 




STANDARD 


STANDARD 

VALUES IN METERS 


RMS DEVIATION 

RMS 

DEVIATION 

CONVENTIONAL 

GPS 

X 

13.0 3.9 

14.2 

2.0 



y 

20.2 3.5 

18.1 

3.5 



z 

6.1 3.0 

7.8 

4.1 

DIFFERENTIAL 

GPS 

X 

4.1 4.1 

3.4 

3.4 



y 

7.8 3.3 

5.7 

4.3 



z 

3.1 1.3 

3.2 

2.8 

Figures 

3-11 

and 

3-12 show the performances for the constant 

case and for 

the adaptive 

case 

, respectively, during the turning segment only. 

Table 3-8 

compares the 

RMS 

and 

standard deviations for both cases. Figure 

3-13 shows 

the process noise history during the turn* 



Table 3-8. 

PVAT Filter Performance (Meters) 

, Adaptive Mode 




Turning Segment Only 






Diagonal Q Matrix 






CONSTANT Q MATRIX 

ADAPTIVE Q MATRIX 




STANDARD 


STANDARD 

VALUES IN METERS 


RMS DEVIATION 

RMS 

DEVIATION 

CONVENTIONAL 

GPS 

X 

9.8 4.2 

14.1 

1.7 



y 

24.0 3.5 

19.9 

4.5 



z 

2.2 0.4 

2.5 

0.7 

DIFFERENTIAL 

GPS 

X 

5.3 4.3 

2.7 

2.4 



y 

11.7 3.7 

7.7 

5.0 



z 

4.8 0.16 

4.3 

0.65 


In this case, the filter seems to track better in the adaptive mode in 
the X direction (smaller standard deviation) but the performance is slightly 
degraded in the other directions. In the differential mode, the RMS accura- 
cies are improved in all axes. 

Overall, though, the improvement of the adaptive scheme is not dramatic 
and is characteristic of the other schemes tried on the GPS program for vari- 
ous receivers. One such method is described in the next section. Later in 
this study, we report on results during landing approach type patterns where 
increased dyn£unics takes place. 

3.3.2 Residuals Tests 


This method was implemented by Magna vox in their Z-set. It consisted of 
calculating four filtered pseudorange and four filtered deltarange residuals 
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and multiplying the diagonal elements in the Cholesky covariance matrix by 
their simple average* Although this method was not repeated, it is similar to 
the one used here in that it increased the covariance (and therefore increases 
the gains) when the residuals get large (indicating poor filter tracking). 
Since the Z-set method is a simpler version with weaker assximptions about how 
the various axes can be combined, the Z-set method will be less responsive 
than the method applied here. 
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Figure 3-1 la. PVAT Filter Perfomance During Turn, X-Axis Constant Q-Matrix 
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Figure 3-11b* PVAT Filter Performance During Turn, Y-Axis Constant Q-Matrix 
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Figure 3-1 1c. PVAT Filter Performance During Turn, Z-Axis Constant Q-Matrix 
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Figure 3-1 2a. PVAT Filter Performance During Turn, X-Axis Adaptive Q-Matrix 
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Figure 3-12b. PVAT Filter Performance During Turn, Y-Axis Adaptive Q-Matrix 
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Figure 3- 12c. PVAT Filter Performance During Turn, Z-Axis Adaptive Q-Matrix 
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Figure 3-13. Process Noise History During Turn 
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IV. MISSIOW-TArU)RKD SATKLLITB SKLECTION 


Much has been written on the subject of ”optimal" satellite selection for 
GPS. Techniques generally optimize the geometry of the four satellites needed 
for stand-alone, continuous GPS tracking. These techniques minimize the posi- 
tion dilution of precision, or PDOP, which is the root sum square of the geo- 
metry-induced errors in each of three orthogonal axes. Minimum PDOP, there- 
fore, results in the minimum sum of the squares of the errors in each axis, 
assuming that the satellite to user range errors are all equal and therefore 
can be normalized out. 

Other techniques have also been proposed, some of which weigh a priori 
known values of each satellite’s measurement error variance [5]. In addition, 
the satellite selection process can weigh observed errors by mathematically 
inferring their source. Although these more sophisticated techniques clearly 
make better use of all available information and will probably be signifi- 
cantly better if the satellites are not uniformly corrupted, they do pose com- 
putational problems. Observation of operational satellite data will decide 
this tradeoff. In any case, geometry-based selection algorithms can be 
studied independently of range variance-based methods, since geometry methods 
simply assume optimal performance of any range variance-based method, and thus 
would augment such a technique. 

Of course, the subject of optimal satellite selection can be avoided 
altogether by employing an "all- in- view” tracking strategy. This technique 
observes all satellites either continuously (with a multichannel or multiplex 
set) or sequentially. There are tradeoffs involved in this case also, how- 
ever, due to the complexity and increased uncertainty of the tracking and 
switching environment created. 

4.1 SATELLITE SELECTION ALGORITHM CONCEPT 

An immediately applicable geometry-based satellite selection concept is 
to consider the mission reo[uirements in choosing satellites. This was inves- 
tigated in the present study. To provide the foundation for this geometry- 
based analysis, however, a brief derivation of the GDOP (PDOP plus the time 
term) concept is first presented. 

GDOP is defined for a system whose measurements, z (the pseudorange 
errors), are related to the error state, x, by the expression: 

z * Hx + V 

where 

V =* unmodeled errors (white noise) 

H = direction cosines to the four satellites 

and 


E[v ^] = R 
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For a best linear unbiased estimate of x given z, 

— T-1 -1 T-1— 

X » [H R H] HR z 

The error in this estimate has the covariance: 

cov (x) = E[(x - x)(x - x)^] * R ^ 

GDOP is defined by assuming that the measurement errors are uncorrelated 
and identically distributed so that: 

R =» 0^1 

Thus, the error covariance is: 

cov (x) = 

where GDOP is defined as: 

GDOP (tr 

Note that the GDOP term is a ” compromise" between the various components of 
the [H^H] ^ matrix, where a less than minimum DOP value in one aocis of a par* 
ticular constellation's GDOP may be sacrificed (accepted) to avoid selection 
of a very large DOP value in some other axis of another constellation* 

The concept to be investigated is whether or not, for a particular mis- 
sion application, one may want to weigh certain axes that are more important 
to mission needs. In the landing situation, the critical axes of concern are 
the vertical and cross-track coordinates. The vertical axis carries the most 
restrictive specification in the landing criteria as shown by the FAA Naviga- 
tion System Accuracy Standards presented in Table 4-1 [6] • The implication 
from the table is that the along-track axis accuracy can be relaxed somewhat. 


Table 4-1. Minimum Guidance Accuracy 


Lateral 


Vertical 


Category 

(ft) 

(m) 

(ft) 

(m) 

(ft) 

(m) 

I 

100 

(30.5) 

30.0 

(9.1) 

10.0 

(3.0) 

II 

50 

(15.3) 

15.0 

(4.6) 

4.5 

(1.4) 

IIIABC 

0 

(0) 

13.5 

(4.1) 

1.8 

(0.5) 

To pursue 

this possibility. 

the GDOP 

matrix was 

modeled 

in a " 


mission coordinate frame" with three orthogonal axes in the along-track, 
cross-track, and vertical directions as shown in Figure 4-1. GDOP, being the 
root sum square of all three coordinates, is the same in either coordinate 
frame, of course. The vertical axis is oriented normal to the local tangent 
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plane, although it could be tilted by the glideslope angle if desired. At a 
typical 3* glideslope angle, the difference would be negligible. The error 
covariance in the landing coordinate frame will be; 

_ -T T -1 -1 

cov(x) - a [A H HA ] 

- a 

where 

2 

* along-track dilution factor 
A 

2 

* cross-track dilution factor 

2 

* vertical deviation dilution factor 

2 

V » time dilution factor 
T 

A * Rotation matrix from normal GPS coordinate frame to A-X-V frame 

A satellite selection algorithm optimized for this fratme may seek to min- 
imize Vy or some weighted combination of Vy and V^. 




A = ALONG-TRACK 
X = CROSS-TRACK 
V » VERTICAL 


Figure 4-1. Landing Approach Oriented Coordinate Frame 
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4.2 SATELLITE SELECTION ALGORITHM SIMULATION RESULTS 


The above GDOP coordinate fraune rotation was modeled in NASA's DIFFGPS 
Simulation OOPS Analysis Module. Then, representative areas of operation for 
the remote helicopter mission were identified, and both conventional minimum 
POOP and the modified satellite selection algorithms were executed over a 
12-hour period. The satellite constellation was selected as the proposed 
18-satellite, 6-orbit configuration [7]. 

The modified satellite selection algorithm used for these runs is an 
even-weighted ”XVDOP”, where the criterion was a sum square of the cross-track 
and vertical error values: 

XVDOP 

In addition, a VDOP criterion is analyzed, which minimizes the error only in 
the vertical direction. 


Figure 4-2 presents a typical plot of dilution of precision over a 
12-hour period. The plot presents the "VDOP" (vertical dilution) component 
values achieved by three different satellite selection criteria. The first 
criterion is PDOP, position dilution of precision, which utilizes all three 
axes and is represented by "XYZ" in the figure. The next criterion is XVDOP 
as described above, represented by Y2 in the figure. The last criterion is 
VDOP, which is the criterion of minimizing only the vertical component, 
regardless of the values in the other two axes (in this case, cross-track 
accuracy may suffer). 


The XVDOP case is dependent on the azimuth orientation of the vertical 
plane, of course. In all cases presented here, the orientation was normal to 
the runway heading of the referenced city. In other analyses, not presented 
here, the sensitivity of satellite selection to azimuth orientation of the 
vertical plane was studied. Although differences did exist, they were gener- 
ally not significant to the VDOP value. Therefore, the random seunple pre- 
sented here is considered "representative." 


DiLirrioM or precision plots for best satellite combos 
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Figure 4-2. Dilution of Precision, Seattle 


The statistical results are calculated as follows. First, the PDOP value 
is calculated. Next, the percentage of time over the 12 hours that VDOP dif- 
fered from the nominal PDOP selection case is calculated. Finally, the per- 
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centage improvement of VDOP, during those periods when it differs from the 
nominal case, is calculated. The results are tabulated in Table 4-2. 


Table 4-2. VDOP Statistics for Seattle, WA Case 


Nominal PDOP 
VDOP 

Percentage of Time Different 
from PDOP Criterion 
Improvement Amount When 

Better than PDOP Criterion 


Selection 

Criterion 

XVDOP 

VDOP 

5.2 

5.2 

22.8% 

32.4% 

.35 

.35 


The improvement in VDOP of .35 represents about a 10% improvement, which 
as indicated, occtirs about 23% of the time for the XVDOP selection criterion 
case and 32% of the time for the VDOP selection criterion case. Table 4-3 
presents results for several other representative locations, for the XVDOP 
criterion only. It should be noted that VDOP is not always better, when dif- 
ferent, in the XVDOP selection criterion case. For these results, it was bet- 
ter on the average over 75% of the time. Better weighting of the cross-track 
and vertical terms would correct this. 

One notable result illustrated by the plot in Figure 4-2 is that this 
improvement occurs primarily when overall DOPs are "good", and not during the 
VDOP "spikes” that occur due to changing geometry. This result is in general 
true for all cases tested; this is an unfortunate result since it would be 
beneficial to find a means to improve VDOP during these "bad" periods. 

A possible complementary solution would be to add an altimeter measure- 
ment and state. The altimeter's inherent inaccuracies would probably elimi- 
nate its influence during periods where the VDOP is "good," but this is where 
the XVDOP selection algorithm improves performance. However, in periods where 
the VDOP is poor and where the XVDOP algorithm was shown to have no effect, 
the altimeter input is likely to be more heavily weighted thereby substituting 
its vertical "measurement" for the poorly resolved GPS vertical observation. 
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Table 4-3. VDOP Results for XVDOP 



Selection Criterion 


Percent 


of Time 

Location 

Different 

Fairbanks 

14.9 

Kodiak 

41.5 

Cold Bay 

19.1 

Juneau 

45.6 

Seattle 

22.8 

Seattle W. 

21.2 

Portland 

58.1 

Portland W# 

70.5 

SFO 

7.9 

SFO W. 

7.9 

L.Ae 

11.2 

L.A. We 

11.6 

Bangor 

18.7 

Bangor £• 

18.7 

St. John 

12.0 

Ste John £• 

12.4 

Average 

17.4 

W. - 100 miles 

west of city 

£• » 100 miles 

east of city 


Improve- 
ment 
in VDOP 

.12 

.24 

.20 

.23 

.35 

.26 

.24 

.17 

.32 

.26 

.23 

.21 

.29 

.25 

.30 

.27 

.25 



V. SENSOR AIDING 


5.1 BAROMETRIC ALTIMETER 


The addition of baroaltimeter inputs provides improved vertical accuracy 
in high VDOP conditions. A PVAC filter was modified to include a baro state, 
and its performance was compared to that of the original PVAC filter. The 
selected baro state is (^Qps “ ^Baro^* 

Two modes are used; the calibrate mode when baro measurements are not 
used (POOP < 5] in the position calculation and the measurement mode, when 
baro measurements are used (PDOP >5)* In the calibrate mode, ^Qp 3 " ^Baro 
estimated by filtering the b^pg - measurements. Tbe baro error equation 

is: 

5h^(k) = (1 - ^ 
b 

with Wjj s N(0, Njj) 

Njj * (9.6 ft)Vsec 
5hj^(0) = 500 ft 
- 10,667 sec 

In the measurement mode, the bg^^^ measurements are processed directly and 
contribute to tbe position solution. 


- 1) + W^(k - DAt 


b« * - (b^^^ - b^ ) 

Baro GPS GPS Baro 


5.1.1 Truth Baro Model 


The baro error model is 


with 


5h*e +e h + cv 
b p hsf sp 
o 


with 


(e ) * (1 - 0 ) At) (e ) + w At 

p k+1 alt P P 

o o o 


w - N(0, Q ) 

P P 

e °(0) = 20 ft 

2 2 
Q = 2u) a = (6.9 ft)/sec 
p alt alt 


-4 -1 

( 1 ) = v/d = 10 sec 

alt alt - 

where 

e = error due to the variation in altitude of a constant pressure 

p 

^o surface 
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6j^g£ * scale factor error due to non-standard temperature 

Cgp * coefficient of static pressure measurement error 

dg 2 ^^ » correlation distance of weather patterns 

=» standard deviation of the variation in altitude of a constant 
pressure surface 

V =« speed 

The baro measurement error consists of a bias with a "small eunount" of 
noise. Let us select: 

R * 25 ft^ 

The process noise is 
with 

Ah = 5 ft 
Av = 1 ft/sec 


or 


Q = (9.6 ft)^ 
5.1.2 Simulations 


The Constant Acceleration Turn pattern used previously was used for the 
baroaltimeter simulations. A forced satellite switch was implemented leading 
to a PDOP history as shown in Figure 5-1. The PVAC filter without baro inputs 
was r\in as a benchmark. Figure 5-2 shows its performance. The baroaltimeter 
filter was run in both the calibrate mode (no baro measurements are processed 
when PDOP is less than five) and in the measurement mode (baro measurements 
are processed throughout). These results are provided in Figures 5-3 and 5-4, 
respectively. There is little difference in performance between the calibrate 
mode and the measurement mode. In the calibrate mode, the baro measurements 
are not used to calculate the position when PDOP is less than five. The baro 
bias error is, however, estimated continuously. The vertical accuracy is 
slightly better in the non-calibrate mode (i.e., when baro measurements are 
always processed in the position determination), due to the additional infor- 
mation in the vertical direction provided by the baro measurements. If these 
measurements are well-modeled in the filter, they will provide increased 
accuracy. It is therefore advantageous to always process them in the calcula- 
tion of the position. 
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Figure 5-1- GDOP History 
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Figure 5-2a* PVAT Filter Performance, X-Axis Degraded POOP 
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Figure 5-2b. PVAT Filter Performance, Y-Axis Degraded POOP 
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Figure 5-2c. PVAT Filter Performance, Z-Axis Degraded POOP 
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Figure 5-3a. ALT Filter Performance, X-Axis Calibrate Mode 
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Figure 5-3b. ALT Filter Performance, Y-Axis Calibrate Mode 
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Figure 5-3c. ALT Filter Performance, Z-Axis Calibrate Mode 
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Figure 5-4a. ALT Filter Performance, X-Axis Measurement Mode 
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Figure 5*4b. ALT Filter Performance/ Y-Axis Measurement Mode 
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Figure 5-4c. ALT Filter Performance, 2-Axis Measurement Mode 
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RMS and standard deviations of the errors in all axes are provided in 
Table 5-1 for both the PVAC filter and the PVAC/baroaltimeter filter (in the 
calibrate mode ) . 


Table 5-1. ALT Filter Performance (Meters), 

Versus PVAC Filter Performance Constant Acceleration Turn Pattern 




RMS 

PVAC 

STANDARD 

DEVIATION 

RMS 

ALT 

STANDARD 

DEVIATION 

CONVENTIONAL GPS 

X 

28e3 

14.0 

13.7 

6.8 


y 

48.2 

11.2 

33.4 

3.8 


z 

377.0 

241.0 

8.9 

7.6 

DIFFEPENTIAL GPS 

X 

KD 

. 

00 

in 

. 

8.3 

7.2 


y 

14.2 

9.3 

4.5 

3.7 


z 

374.0 

237.0 

8.4 

7.2 


Clearly, the performance of the ALT filter is superior, especially in the ver- 
tical direction. Differential GPS helps remove the common biases in the mea- 
surements for both filters. It is not improving much vertical accuracy of the 
PVAC filter. This is because the common biases in this direction are small 
and the geometry is poor (large VDOP). Range errors will therefore translate 
into large altitude errors. For the ALT filter, the altitude solution will be 
heavily based on baromeasurements and will therefore be much more accurate. 

In the PVAC performance to be better, the process noise should be decreased 
when the OOP is increased so as to rely more heavily on a priori Information 
and less on new GPS measurements. This would be more in line with the correc- 
tor filter tuning and should provide for better error cancellation in the dif- 
ferential mode. It is clear, however, that the inclusion of baroaltimeter 
inputs dramatically improve performance without the need for adaptive process 
noise features. 

The ALT filter performance is studied further in the next section using a 
landing type profile. 

5.2 INERTIAL COMPONENTS 

Although inertial sensors were not modeled and simulated in this effort, 
because of their potentially large contribution to an integrated GPS landing 
guidance system, conceptual analysis was performed. Use of inertial sensors 
could substantially change the modeling used by the Kalman filter processing 
the GPS measurements, for example. 

The inertial sensors considered here are of the "low cost" variety, that 
is, not including complete stabilized platform (inertial navigation system) 
configurations. Instead, the emphasis is on isolated accelerometers and 
gyros. 
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Since the major objective in a GPS-based landing guidance system is to 
reduce the vertical component of error/ the logical choice is to incorporate a 
vertical accelerometer. The vertical accelerometer outputs precise changes in 
velocity/ so it would provide an excellent source of information about motion 
in a "vertical** direction. Roll and pitch motions of the aircraft would dis- 
tort this output as a measure of local level vertical/ but since only the 
landing domain is of interest in this design/ such effects would be minimal. 
For example/ a 10* bank angle (reasonably severe maneuver in a typical final 
approach) would resolve a .2 g local level axis-vertical maneuver (again 
reasonably severe) into an erroneous .203 g vertical maneuver. Similarly/ a 
flat 10* bank turn maneuver would be interpreted as a little over 1/100 of a g 
climb. Appropriate modeling and weighing in the Kalman filter would not 
respond to these small values for short durations. 

To incorporate a single vertical accelerometer in an integrated GPS Kal- 
man filter/ several options are available. The accelerometer error states 
(e.g./ bias and scale factor) could be modeled in the state vector/ enhancing 
the accuracy of the delta velocity information while increasing the dimension- 
ality of the processing and therefore processing load. Such modeling is fre- 
quently used in full INS/GPS integration/ but the major objectives for such 
applications are long term stability/ constant calibration for GPS outages/ 
and stability of information under maneuvering conditions for precise carrier 
loop aiding (doppler prediction). For this application/ most of these objec- 
tives are not relevant. 

Alternatively/ the small accelerometer errors which impact this applica- 
tion only in a short term sense can be ignored and handled by a small amount 
of process noise. The delta velocity output is used to propagate the vertical 
velocity state of the filter between measurement updates instead of the accel- 
eration modeling studied in the earlier sections of this report. Note the 
significance of this concept; the acceleration measurement from the accelero- 
meter is likely to be far superior to any a priori model or computation/ and 
yet the state vector can actually be reduced in size by eliminating the verti- 
cal acceleration state (since it is known so need not be estimated). In fact/ 
it may also be inconsequential to model lateral accelerations as well given 
that the vertical axis of concern is well modeled. 

Extending this concept/ a couple or triad of accelerometers could be used 
to measure the other axis or axes. Small mounting errors could be adequately 
handled by a small amount of process noise. 

Use of other on-board instruments could be made to further resolve the 
accelerometer outputs/ especially if performance in dynamics is important. 

The attitude indicator provides a fairly good steady-state bank and pitch 
attitude. This would be easy to integrate if digital outputs were available 
with sufficient resolution/ latency/ and time tag accuracy. 

These various concepts provide several possible low cost enhancements to 
stand-alone GPS operation. Their relative benefits and practical significance 
must be determined by simulation and eventual flight test. 
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VI* APPROACH AND LANDING SIMOIATIOW 


To create a realistic landing scenario for the DIFFGPS Monte Carlo Simu- 
lation, a landing path was devised that emulated an actual manual helicopter 
landing approach. The objective was to accurately model the dynamics of a 
manual approach. Accuracy of the dynamics was important since the GPS filter 
modifications were designed to model essentially constant acceleration dyna- 
mics which may or may not be true of the actual approach path. 

6.1 DEVELOPMENT OF SIMULATED LANDING PATH 

An accurate simulation could be produced by explicitly modeling a heli- 
copter, instrument landing system, and closed loop pilot, of course. However, 
it was felt that the complexity of such a model, along with the difficulty in 
tuning the aircraft and pilot models, was excessive for this effort. 

Instead, an auialysis was made of manual helicopter landing data and a 
mathematical description of the profile was devised. Figure 6-1 shows a sam- 
ple of the landing data from Reference 5. These tracking data were gathered 
during NASA Beacon Landing System testa on an SH-3 helicopter. Other example 
data included helicopter MLS Approaches and ctirved approach testa. Analysis 
of the data in Figure 6-1 reveals a basic oscillatory pattern about the glide- 
path, particularly in the lateral aucis. The vertical axis exhibits some of 
these oscillations although less pronounced and more heavily damped. 
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Figure 6-1. Actual Helicopter Landing Tracking Data 
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Based on these data/ a ** randomized** sinusoidal path motion was selected 
in the lateral and vertical axes* The form of the model (for the lateral 
case) was: 

y(x) - Ae sxn(Cxe + E) + F 

where 

A, B, C, 0 are exponentially correlated random variables 

E is a constant randomly initialized 

F is a Gaussian random variable for gust term 

Note that both the period and the eunplitude of the sinusoid were allowed 
to randomly adjust, of course ensuring that the path was nevertheless continu- 
ous* The period and amplitude of the motion were also designed to generally 
reduce in magnitude as the approach neared touchdown, again consistent with 
the appearance of the data* 

To calculate the parameters in the models, coarse measurements were made 
of the period and amplitude of the glidepath deviations at each end of the 
approach for both the horizontal and vertical axes* These data are shown in 
Table 6-1* Note that any bias in the glidepath deviation was ignored; it is 
not significant in determining aircraft dynaunics during approach for small 
values* The final lateral and longitudinal path equations were: 

f(x) * Ae®^ sin(Cxe^^ + E) + F 


For A, B, C, D, E: 


X » xe 


-G/G A+ 
o 


+ a X 


where the values are: 


a 

-A 



-26/G A+ 
e ° N(0,1) 


Lateral 


80 

90 

1 

3040 


Vertical 


17 

21.4 

1 


3040 


B .0000153 .0000228 

0_ .0000031 .0000047 

0 
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Lateral 


Vertical 



1 

1 

Sb 

50 

50 

c 

.0032 

.004 

0 

c 

.00064 

.0008 


1 

1 

®oc 

1500 

1500 

5 

-.0000224 

-.0000224 

‘^o 

.00002 

.00002 

«o 

1 

1 

®0D 

3040 

3040 

B 

u(-3.14,3.14) 

u(-3.14,3.14) 

F 

N(0,1) 

N(0,2) 


Velocity was also modulated along the path to approximate a pllot/alr- 
craft response when maneuvering vertically. This was modeled as "speed along 
the path” for trajectory propagation purposes: 

V(x-1)+N(0,.5) X > .750 

V(x) - V(x-1)+32.2 ^ •‘^50 

\ X* 1 / "X 

V « 135 ft/sec 
o 


where : 


X Is distance along path 

D Is total approach path distance (i.e., glideslope starts at .750) 
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While it can be argued that this model is not an exact representation o£ 
the observed data/ the fact that the match is approximate/ the dyn<unics have 
infinite derivatives/ and the randomness induced should provide more than ade- 
quate "realism" for the Kalman filter attempting to follow this motion. 

Figures 6-2 and 6-3 present plots of 10 approaches in each axis which can be 
favorably compared with the recorded actual data of Figure 6-1 (again noting 
the absence of glidepath bias errors in the simulated data). 


Table 6-1. Path Dynamics for Landing Data 

Beginning of Approach End of Approach 



Peak- 

Peak Oscil- 
lation 

Standard 

Deviation 

(2s) 

Period 

Peak- 

Peak Oscil- 
ation 

Standard 

Deviation 

(2s) 

Period 

Vertical 

102 

156 

3541 

47 

63 

1672 

Horizontal 

313 

406 

4459 

188 

141 

2108 
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Figure 6*3 • Horizontal Axis Simulated Landing Approaches 
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6.2 LANDING SIMULATION FILTER PERFORMANCE 


The flight path described in the previous section was used to test the 
filters designed earlier in this study, namely, the PV, PVAC, PVAT, and ALT 
filters. The filters were tuned for a level Constant Acceleration Turn pat- 
tern and their performance was descussed in Sections 3 and 5. In the landing 
simulation, the vehicle experiences more dynamics (of a sinusoidal nature) in 
all axes. Although we found out that the PVAC filter was probably the best 
choice for the Constant Acceleration Turn pattern, it is possible that the PV 
filter will perform just as well for the landing path. This is because of the 
continuously changing acceleration (in a sinusoidal way). Modeling of the 
acceleration as constant or even centripetal will result in filter lagging and 
could degrade performance as compared to the PV filter. In the next sections, 
performance of the filters is compared, and the results of Monte Carlo runs 
are provided. 


6.2.1 PV, PVAC, PVAT Filters Comparisons 


The landing profile described in Section 6.1 was flown with the PV, PVAC, 
and PVAT filters. Three cases were tested for the PVAT filter; constant pro- 
cess noise matrix Q, variable matrix Q, and adaptive matrix Q. Results are 
provided in Tables 6-2 through 6-4 and in Figures 6-4 through 6-8. 


The PVAC filter, as for the constant acceleration turn, out performs the 
PV filter in both the conventional and differential modes. The advantage is 
most pronounced in the East direction (lateral motion) with a 44% improvement 
in standard deviation in the conventional mode and a 35% improvement in the 
differential mode. At the RMS level, a 30% improvement is achieved in the 
East direction for the differential mode, with only a 5% improvement in the 
conventional mode. In the North directions, a small improvement is achieved 
with the PVAC while in the vertical direction the performance is even. 


The various PVAT filters perform pretty much identically. Their accuracy 
is similar to that achieved with the PVAC filter. Although the PVAC 
performance is better than that achieved with the PV filter, the improvement 
is not as large as was obtained on the Constant Acceleration Turn pattern. 

This is because the acceleration vector is continuously changing (in a 
sinusoidal way) for the landing pattern, and the estimate of acceleration will 
therefore always be lagging, thereby reducing the performance improvement over 
the PV filter. 


Table 6-2. PV Filter Performance (Meters) Landing Profile 

RMS Standard Deviation 



East 

17.3 

6.6 

Conventional GPS 

North 

17.3 

6.8 


Up 

7.1 

4.0 


East 

7.5 

6.5 

Differential GPS 

North 

9.5 

7.8 


Up 

4.2 

4.2 
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Figure 6-4. PV Filter Performance Landing Profile 
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Table 6-3. PVAC Filter Performance (Meters) Landing Profile 




RMS 

Standard Deviation 


East 

16.5 

3.7 

Conventional GPS 

North 

16.6 

6.5 


Dp 

7.1 

4.1 


East 

5.2 

4.2 

Differential GPS 

North 

8.7 

7.3 


Up 

4.2 

4.1 

Table 6-4. 

PVAT Filter 

Performance 

(Meters) Landing Profile 


RMS Standard Deviation 




Constant Variable Adaptive 

Constant 

Variable Adaptive 



1 

£ 

£ 

£ 

£ 

£ 


East 

16.5 

16.5 

16.4 

3.7 

3.7 

3.5 

Conventional GPS 

North 

16.8 

16.9 

16.7 

6.7 

6.5 

6.5 


Up 

7.1 

7.1 

7.7 

4.1 

4.1 

4.6 


East 

S.2 

5.2 

5.1 

4.2 

4.2 

4.1 

Differential GPS 

North 

8.8 

8.8 

8.5 

7.4 

7.2 

7.0 


Op 

4.2 

4.2 

4.5 

4.1 

4.1 

4.4 


In order to generate more meaningful statistical data on the filters per- 
formance, ten Monte Carlo runs were performed on the PV, PVAC, and PVAT (Con- 
stant Q) filters. Results are reported in Section 6.2.3. 
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Figure 6-6. PVAT (Constant Q) Filter Performance Landing Profile 
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Figure 6-7. PVAT (Variable Q) Filter Performance Landing Profile 
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Figure 6-8. PVAT (Adaptive Q) Filter Performance Landing Profile 


43 - 3020-83 




6* 2. 2 ALT Filter Performance 


The ALT filter developed for the Constant Acceleration Turn pattern was 
run over the landing profile and for the same degraded DOP condition as was 
used for the Constant Acceleration Turn pattern. As a comparison, the PVAC 
filter was run under the seune conditions* Position errors are shown in Figure 
6-9 for the PVAC filter, Figxire 6-10 for the ALT filter, and Figure 6-11 for 
the reference receiver filter. The GDOP history is shown in Figure 6-12. 
Performance comparisons are shown in Tables 6-5 and 6-6* 

The performance degradation at the GDOP transition is clearly visible for 
the PVAC filter as the errors grow very large, especially in the vertical dir- 
ection. The degradation is more subtle for the ALT filter, and although the 
altimeter helps substantially at the GDOP transition, the vertical accuracy 
degrades to 100 m error. Additional filter tuning in the vertical direction 
is probably required to enhance further the bauro inputs over the GPS measure- 
ments and the system modeling. 
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Figure 6-9. PVAC Filter Performance, Degraded DOP, Landing Profile 
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Figure 6-10. ALT Filter Performance, Degraded DOP, Landing Profile 
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Figure 6-11. Reference Receiver Filter Performance, Degraded OOP 
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Figure 6-12. Landing Profile and GDOP History 
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Table 6-5. PVAC Filter Performance (Meters) Landing Profile 

Degraded OOP 




RMS 

Standard Deviation 


East: 

18.3 

12.0 

Conventional GPS 

North 

41.3 

11.9 


Up 

220.4 

212.7 


East 

12.4 

9.6 

Differential GPS 

North 

8.2 

7.5 


Op 

218.1 

167.0 


Table 6-6. 

ALT Filter 

Performance (Meters) 

Landing Profile 



Degraded DOP 





RMS 

Standard Deviation 


East 

15.3 


9.1 

Conventional GPS 

North 

37.1 


4.7 


Op 

50.9 


38.1 


East 

6.6 


5.0 

Differential GPS 

North 

3.3 


3.1 


Op 

53.9 


36.9 


6.2.3 Monte Carlo Runs 

Ten runs were performed for the PV, PVAC, and PVAT (Constant Q) filters 
using ten randomly generated seeds to drive the error models. Statistics were 
generated for each filter as follows. A mean and standard deviation were cal- 
culated for each run. An overall mean was then generated as the mean of the 
Individual mean values. An overall stauidard deviation generated as the square 
root of the mean of the sq[uares of the individual standard deviations. An 
overall SMS value was generated as the RSS of the overall mean and standard 
deviation. Results are listed in Tedsles 6-7 through 6-9. 

In the conventional mode, the PVAC filter out performs the PV filter in 
both RMS ^md standard deviation, with a 45% improvement in East st 2 uid 2 urd devi- 
ation (which pretty much corresponds to lateral motion) and a 12% Improvement 
in North standard deviation. There is a 1% degradation in vertical standard 
deviation. The RMS values show smaller improvements, 22% in the East, 1% in 
the North, 1% in the vertical. 

In the differential mode, performance comparisons are similar, with the 
PVAC out performing the PV filter in standard deviation by 44% in the East and 
8% in the North. There is a 3% degradation in the vertical direction. For 
the RMS values, the PVAC has an advantage of 42% in the East auid 9% in the 
North, but there is a 3% degradation in the vertical direction. 
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Performances of the PVAC and PVAT filters are virtually identical. The 
results above indicate that it is advantageous to estimate the acceleration in 
the filter* However, the model itself is of little importance, with a simple 
constant acceleration scheme performing just as well as more sophisticated 


models* In the vertical 

direction, however, there 

was 

actually a small degra- 

dation of accuracy when acceleration was estimated. 

The vastly improved lat- 

eral performance should. 

however, make the PVAC filter 

the preferred choice 

for the landing application. 



Table 6-7. PV 

Filter Performance (Meters) Landing Profile 


10 Runs Statistics 




RMS 

Standard Deviation 

X 

9.0 


6.8 

Conventional GPS y 

7.3 


5.7 

z 

31.9 


6.7 

X 

6.9 


6*8 

Differential GPS y 

6.6 


6*5 

z 

7.S 


7.5 


Table 6-8. 

PVAC 

Filter 

Performance (Meters) Lamding Profile 



10 

Runs 

Statistics 




RMS 

Stemdard Deviation 


X 


7.0 

3.7 

Conventional GPS 

y 


7.2 

5.0 


z 


31.7 

6.8 


X 


4.0 

3.8 

Differential GPS 

y 


6.0 

6.0 


z 


7.7 

7.7 


Table 6-9. 

PVAT 

Filter Performance (Meters) 

Landing Profile 




(Constant Q Matrix) 





10 Runs Statistics 





RMS Standard Deviation 



X 

7.0 

3.7 


Conventional GPS 

y 

7.1 

o 

• 

in 



z 

31.7 

6.8 



X 

3.9 

3.3 


Differential GPS 

y 

6.0 

6.0 



z 

7.7 

7.7 



c 
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VII. IANDIWG GLIPgPATH STEERING GUIDAMCS 


This part of the effort developed navigation steering guidance for a 
helicopter executing a repeating touch-and-go landing pattern, including 
scaled lateral and vertical steering guidance during final straight-in 
approach. 

7.1 PATTERN DESIGN 

The landing pattern, shown in Figure 7-1, is designed to meet NASA flight 
test procedures. Settable parameters are listed in Table 7.1. Horizontal and 
vertical steering gains for the non-approach segments of the profile are set 
for nominal path guidance* 

The approach gains were designed to maximize compatibility with typical 
pilot procedures. That is, the approach gain schedule is a function of range 
to glidepath intercept point (GPIP). The maximum extent of the display sensi- 
tive area converges as the aircraft approaches the GPIP in both the vertical 
and horizontal axes, much like the traditional glideslope and localizer sig- 
nals of the Instrument Landing System (ILS). With GPS, of course, there is no 
reason that the "beamwidth” must converge down the glideslope; a constant gain 
"tube" could be implemented. 

However, there are inherent advantages in a converging sensitivity bound- 
ary. In the early stages of the approach, the pilot is capturing the path in 
both the horizontal and vertical auces. He is establishing stable path track- 
ing. Finally, since the aircraft is still high above the ground, the pilot is 
usually less precise in his tracking and unvoncerned about small path devia- 
tions. Therefore, a wide sensitivity area, or low display gain, is appropri- 
ate* 


At the lower end of the approach, the situation is different, particular- 
ly in low visibility conditions. Here it is critical to maintain ground 
clearance and stable approach conditions, and tight path tracking is essen- 
tial. Increasing the sensitivity gives more precise path deviation informa- 
tion to the pilot, who is now controlling the aircraft beam deviations to 
finer tolerances. 

The other major feat\ire of this design is separate specification of cap- 
ture, mid-range, and final path sensitivity ( "beamwidth" } boundaries. This 
enadDles low gain on the display while capttiring the horizontal and vertical 
paths, and higher gain once stable path tracking is established. Furthermore, 
the various sensitivities are selectable in flight so that successive 
approaches can optimize the gain selections. 
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Table 7-1. Settable Landing Pattern Parameters 


Parameter 

Description 

Default 

PA 

Pattern Altitude 

610 ni • AGL 

DL 

Downwind Length 

8 nm 

BL 

Base Length 

2 nm 

AFD 

Approach Final Distance 

5 nm 

XCV 

Max Lateral Deflection 

• 5 nm 

YCV 

Max Vertical Deflection 

152.4 m 

ETR 

Expected Turn Radius 

• 5 nm 

CR 

Capture Range 

2.5 nm 

CLW 

Lateral Capture Window 

1219.2 m 

XMLW 

Middle Lateral Window 

61.0 m 

GPIPLW 

GPIP Lateral Window 

3e0 m 

CVW 

Vertical Capture Window 

61«0 m 

XMVW 

Middle Vertical Window 

15.2 m 

GPIPVW 

GPIP Vertical Window 

6*1 m 

GLIDE 

Glideslope Angle 

6 deg 

APBEAR 

Approach Bearing 

11 deg 

VF 

Velocity Factor 

100 

SF 

Scale Factor 

2.5 

GLAT 

GPIP Location/ Crows Landing 37* 

24* 47.79" 

GLON 

GPIP Location, Crows Landing- 121* 

06’ 30.39* 

GALT 

GPIP Location, Crows Landing 

5.4 m 


(Ellipsoid 

Height) 


Note: 


1 • XMLH Location • XMVW Location 
2. Middle Window Dist - GPIP « AFD-CR 
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Figure 7-1. Landing Pattern 
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7.2 APPROACH LATERAL STEERING ALGORITHM 


The lateral steering algorithm was formulated as basically a proportional 
plus derivative control strategy. This formulation assures that pilot steer- 
ing commands are always "toward" the landing direction. Closed commeuids by 
the pilot results in "exponential" capture of the path (solution to the homo- 
geneous first order differential equation). 

A "design" approach speed of 60 knots was selected. Although this para- 
meter could have been specified as a variable using actual total helicopter 
velocity, it was decided that the relationship of aircraft control and iner- 
tial path tracking was better maintained by using a design velocity. Thus 
instead of always comm£mding the same path in space regardless of velocity, 
the algorithm will consnand later emd sharper turns at lower velocities. 

The gain strategy was designed to handle ILS unit in use on the test air- 
craft. Scaling was 0 to 250, fly left, 250 to 500, fly right; 0 to 250 fly 
up, 250 to 500 fly down; and 0 to 200:0 to 20.0 nautical miles rauige to GPIP. 

The derivation of the lateral steering algorithm follows: 



250 , 

Uy 


Form: 

C » k^(y + Y y) * ^2 
C » Steering command 
Boundary conditions: 

A. When y « y^ax' y * 0 

Then C » 0 (Full scale left turn) 

B. When y ■ y = 0 

Then C <■ 250 (Centered display) 

C. When y " 0, y « -•■30.48 m/s 

Then C « 0 (Full scale left turn) 


NOTE: Design approach speed 


60 kts (30.89 m/s) 
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From (A) ; 


° ^max^ " ^2 

From ( B ) : 

250 - ^2 
From ( C) : 

0 » le^(30.48) + J«2 
Combining (2) and (3): 


250 


- 8.2 


1 30.48 

Combining ( 1 ) , ( 2 ) and ( 4 ) 

0 - -8.2 4 + 250 

T nuuc 

T * .033 y 

•^max 

Final Expression: 

C » -8.2 [y + + 250 

l^maxl 


( 1 ) 


(2) 


(3) 


(4) 


(5) 


( 6 ) 


where the absolute value sign handles the general left or right of path case. 
2nd Condition: 

For X ■ 0, y < 0 -> 250 <C< 500 
For x*0,y>0->0 <C< 250 


NOTE: y,^,, = f(x) 


The algorithm is propagated over 2.4 seconds at 0.1 second intervals. 
Normal updates of the 2.4 second propagation occur at the GPS measurement 
update rate, typically, 1.2 seconds. 

7.3 APPROACH VERTICAL STEERING ALGORITHM 

The approach vertical steering algorithm is identical to the approach 
lateral steering algorithm with three major exceptions. The gains are struc- 
tured to correspond to the vertical axis parameters (YCV, CVW, XMVW, AND 
GPIPVW), the scaling corresponds to the Instrument vertical scales, and the 
vertical velocity damping is adjusted for a non-zero nominal vertical velocity 
(due to the descent). 
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The derivation of the vertical steering algoritha is as follows: 
From: 

C-)c^(z' + ^z’) +^2 
C ■ Steering command 
Boundary Conditions: 

A. When z* ■ “ 0 

Then C • 500 (Full scale fly down) 

B. When z' ■ -z' “0 

Then C * 250 (Centered Display) 

C. When z' ■ 0/ z'» -6.4 m/s 

Then C « 500 (Full scale fly down) 

MOTE: Design approach speed 60 knots (30.89 m/s) 

z* « z - 2(3 i,ide 

2 2 

z' “ z - tan (Glideslope) x + z 


From (A) : 


--N^Y *’«ax^ *’'2 


From ( B ) : 


250 


From (C): 


500 - -lc^(6.4) + yt^ 

Combining ( 2 ) euid ( 3 ) : 

250 

’'i 

Combining ( 1) « { 2 ) , and ( 4 ) 

500 » 39.1 (^ z’ ) + 250 
Y max 


.156 z' 


max 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 
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Final Expression: 


C » 39.1 [2* + (-7^) 2'] + 250 (6) 

I ^ max I 

where the absolute value sign handles the general up or down of path case. 
2nd Condition: 

For X » 0, 2 '< 0 -> 0 <C< 250 

For X » 0, 2 ' > 0 -> 250 <C< 500 

NOTE; 2 * » f(x) 

max 


The algorithm is propagated over 2 . 4 seconds at 0 . 1 second intervals . 
Normal updates of the 2.4 second propagation occur at the GPS measurement 
update rate, typically, 1.2 seconds. 


7.4 PATTEPN STEERING ALGORITHMS 

The steering algorithm for the "non-approach" segments of the pattern, 
namely the ascent, downwind, and base legs, are constant gain horisontal auid 
vertical steering rules. 

The nominal horirontal steering algorithm is based on a maximum instru- 
ment deflection at ±0.5 nm (Max Lateral Deflection). The steering algorithm 
is of the same form as the approach lateral steering algorithm: 

C » -8.2 (y + (— •— ■) y] + 250 
^max 

C s Lateral steering command 
^max * 

In the pattern algorithm case# of course ^ is constant. 

For turn to the next leg, the lateral algorithm switches to steering for 
the next segment at a distance of 0.5 nm (Expected Turn Radius) from the "to** 
waypoint. Distance is calculated from only the x (along track) component of 
the total distance. 

For turn to localizer, the same event occurs except at the intersection 
of the max lateral distance and the aircraft's x distance from the comer way- 
point within the lateral capture area. This turn distance will typically be 
greater than the nominal pattern value of 0.5 nm. 
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In the vertical aucis, the pattern steering algorithm is also of the same 
form as the vertical landing algorithm^ except that vertical pattern maneuver- 
ing is symmetric about z*0 instead of z**0. Specifically/ the pattern steer- 
ing algorithm is: 

C » 39.1 [z + z] + 250 

z 

max 

C * Vertical steering conmand 


*max 


152.4 m. 


7.5 SIMULATION EVALUATION 


Figure 7-2 shows a sample landing approach path in the runway coordinate 
frame. Figures 7-3 and 7-4 show the steering comm^md computation. Figure 7-5 
is em enlargement of the turn portion of the command. 
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Figure 7-2. Path of Simulated Landing Approach 
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Figure 7-3. Calculated Steering Commands 
(Runway Coordinate Frame, Scaled to ILS u Amps Inputs) 

First 275 Seconds 
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Figure 7-4. Calculated Steering Commands - 225 Seconds to Touchdown 
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Figure 7-5. Expanded View of Calculated Steering Commands 
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VIII* COWCLnSIONS 


Several areas were investigated for their potential contribution to 
improving vertical accuracy for a rotorcraft using differential GPS during a 
landing approach. Continuous deltaranging was studied and the potential 
improvement achieved by estimating acceleration was studied by comparing the 
performance of several filters; a position- velocity (PV) filter, a position- 
velocity-constant acceleration (PVAC) filter, and a position- velocity- turning 
acceleration (PVAT) filter. The filters were tuned for a Constant Accelera- 
tion Turn, horizontal profile, then tested on a "rough" landing profile. 

The performance on the Constant Acceleration Turn pattern was best for 
the PVAT filter, but the improvement over the PVAC filter was negligible. The 
PV filter performance was inferior (factor 2 or 3 in standard deviation). The 
same conclusions, although to a lower extent, were also reached for the land- 
ing profile. Several different implementations of the PVAT filter were 
studied. They consisted of a constant turn rate (and process noise matrix), a 
variable turn rate (and process noise matrix) and an adaptive process noise 
scheme. Very minor improvements, not justifying the added processing burden, 
were achieved with the more complex implementations. 

Overall statistics were generated for ten Monte-Carlo runs of the PV, 
PVAC, and PVAT (constant) filters. Results showed the PVAC filter to be the 
most efficient with the PVAT filter performing equally well, but being more 
complex. Vertical performance, however, was not significantly different (1 to 
3% differences) among the various filters. The largest improvement was along 
the cross-track direction. Satellite selection algorithms based on vertical 
errors only (VDOP) and even-weighted cross-track and vertical errors (XVDOP) 
were tested. The resulting VDOP, when using the above criteria instead of the 
more common PDOP criterion, were different 32.4% of the time (VDOP) criterion) 
and 22.8% (XVDOP criterion). Improvement was achieved with the XVDOP criter- 
ion 75% of the time. One problem with the XVDOP criterion is that improvement 
was generally achieved during good DOP conditions, but not during degraded DOP 
conditions. By adding an altimeter, vertical accuracy should be improved d\ir- 
ing poor DOP conditions, and the XVDOP selection criterion would provide good 
accuracy during "good" DOP conditions. 

The inclusion of an altimeter was studied by modifying the PVAC filter to 
include a baro bias estimate. Tests were r\in both on the Constant Accelera- 
tion Turn pattern and the landing approach and clearly showed the improved 
vertical accuracy during degraded DOP conditions, although additional tuning 
in the vertical direction will be necessary to acquire the required level of 
accuracy during landing. 

Flight Test results for raw differential results, not including the 
effects of filter improvements studied in this simulation effort, indicated 
that the differential performance significantly improved overall navigation 
accuracy. Interestingly, the residual errors were non-dominated by geometry- 
independent errors such that GDOP effects were no longer obvious in the 
remaining error plots. 

A landing glidepath steering algorithm was devised which exploits the 
flexibility of GPS in determining precise relative position. A method for 
propagating the steering command over the GPS update interval was successfully 
implemented. 
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APPENDIX A 


and Q Matrix Derivations 


Below are derivations for the and Q matrices for the various filters 
used in OIFFGPS. 

Q^/ Qj,/ Qf represent the power spectral density matrices (or scalars) 
for the noise in the 3 velocity equations, 3 acceleration equations, clock 
phase equations, and clock frequency equations, respectively. 

The is^lementation consists of having the program user select thepower 
spectral densities (or their square roots) for all equations, and letting sub- 
routines calculate the process noise matrices for the discrete implementation. 

Note that these matrices are constant, except for the turn dynamics fil- 
ter, where it is a function of <i>. 

1.0 <(> MATRIX DERIVATIONS 

1.1 NO ACCELERATION STATES 

X ■ V + noise 

6 » f + noise 

V « noise 

f “ -1/T f + noise 

F matrix is 


X V b,f 
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n2 CONSTANT ACCELEKATION 

The cloclc equations are the same as for Section 1.1 and are left out 
here. The equations are: 

X =■ V 

V - a 

a « o + noise 

F matrix is 



♦ is L“^(SI-F)"^ 
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1 . 3 TUBNING SEGMENT 


The clock equations are the same as for Section 1 . 1 and are left our 
here. The equations are; 

X * V 

V * a 

a * -(D^v + noise 

or F matrix is 
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(SI - F)' 


( ) 


2 2 
(s 



2 2 * 
a I 


2 2 
s 


0(At) 


sin^t . I r^t sin<*iAt ^ 

ST“ ^ 1 S“ ^ 


cosuAtl 


sinOJAt 


0 ! -usinuAtl ! cosidAtl 


d(Ai:) 

o “ 


^ cos<»At 

U.2 o 


(cos<»>At - 1) 


For (» very small/ 0 becomes: 


0 I 


which/ as expected/ is the constant acceleration matrix. 
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0 matrix is 


^ j 

1 

— 1 

1 

I , 

sinoiAt ^ 

1-cos(ijAt j 

1 

° 1 

01 ' 
1 

2 

0) 

cosoiAt I 1 

sinojAt ^ 

1 

1 

U) 

0 1 

-o)Sinu)At I ^ 

cosoiAt I 


V X a 
01 » - 
Iv|2 



Derivation of cj 

Assume central acceleration only. Coordinated turn. 



and. 


o)/a 


U) 

va 


a * ojv 



(J 



2 

V 


V X a 

2 

V 
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2.0 Q MATRIX DERIVATIONS 


2.1 NO ACCELERATION STATES 

The matrix was calculated as 


♦ - 



O 


O 


1 tfi 

0 

T 


Q matrix is 


0 

1 

1 

_ U - 

1 

__°_L 

1 

0 

0 

1 

1 

Qvi 1 

0 

0 

1 

1 

1 

- - 

0 1 
1 

^2^ ~0 
0 Qf 


Qj, - Qo’^dt 


.At 


Q^At I 


Qvl 


^ ~ - U 


Qj, QfAt 
0 


O^dt 


.At 


Q At“ I 

.1 — ^ 


Q At 
*v 


Q^At I O 


_£v j ° I 


QfAt^+QjjQ^Ate"^^/"^ 


QjAt^-At/t QfS 


-2At/t 
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p 1 

[ 

3 ^ 

Q At-* , 

V 1 

2 At^ 

‘ n 

__L_4- 

I 

^ L 

1 

1 

_ _ 

1 

° 

1 

1 

0 

1 

L 1 

^ 1 n At 0 

0 1 3 ;V'9f . 

I 2/ r 


The O^s below are the square roots of the noise P*S.D. for the applicable 
equations (i«e«, for the x velocity equation). 


2(1) 

V 3 

2(2) - 2(4) - 2(5) - 0 

Q(3) 

X 

-a 2^tV" 

V 

2(7) - 2(8) - 2(9) - 0 

2(6) 

y 

.0 2^tV 

V 

Q(11) - 2 AtV ^ 


z 

X 

QdO) 

- 2 2 
-»b - 

2(12) » 2(13) * 2(14) » 0 


Q(15) - ^ At 0 - Q(16) » Q(18) - Q(19) - Q(20) 

X 


2(21) 

■ a ^ At 

V 

2(17) 

-0 2,^2 2 
V 


y 


y 

2(28) 

■ a ^ At 

V 

z 

2(22) 

- 2(23) - 2(25) « 2(26) - 0 - 2(27) 

2(36) 

• At 

2(24) 

- a 2 i^2^ 2 

V 




z 



2(29) 

- 2(30) - 2(31) - 0 


0 - 2(33) - 2(34) - 2(35) 

2(32) - 0^^ (At^) 
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2.2 CONSTANT ACCELERATION 

The clock equations are the same as for 2.1 and are not repeated here. 
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«k 



0 


0 


I 

1 

1 

1 


0 


0 


Q At^ 
a 

2 


Q At 


I 

Atl 




2a 




I 


1 

«k • “^0 

Q~ 1 

“a 4 1 

1 

^a 2 

1 

1 “a 2 

1 


sM 1 
2 1 

Q At^ 
a 

1 

8.4t 

Q 

* J 


_«.¥ i 

Q^At 

Sk - 

r At5 1 

20 i 

1 

o 

®a 8 

j 

1 “ 
1 , 

1 Q^AtVe 

1 


' 

^a 8 1 

At^ 1 

®a 3 1 

1 , 

1 Q^At^/2 


i 

'^a 6 1 

«a¥ ' 

a 2 1 

1 

1 s^it 
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2 . 3 TUBNING DYNAMICS 


The clock equations are the same as for 2.1 and are not repeated here. 
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~ - lainoJAt 

,1 

Qj^(1,2) - -jjj I- (l-cosojAt) - 
®a , 3 

= -TT ^ -r - COSOJAt + 
iiT 4 

2jc^^'3) - {sinuiAt - - 

2,^(2, 2) - ^ - -J sin20>At} 

Qa 

2^(2, 3) - ^ (1-cos2uAt) 

2 ^ 

- -- 3 " (1-co82uAt) 

4 w 

“IT i sin2toAt} 

0 ) << £ 

If then Qjj is same as for constant acceleration case. 

/o^sin2(u)At) d(At) * /o^^sin2((i)At) d(uAt) 

- s 


sin2 (ridt 1 
* 4 ^ 

4— ( 1-cos2(oAt) } 
4u) ‘ 

•7 cos2ci)At} 

4 

•“ sin2uAt} 


/q^cos 2 (uAt) d(At) = — cos2(o)At) d(uAt) ■ [~^ + T sin2ci)At] 
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/q^( 1-coswAt)^ dAt » *2 

-0)At . 4 d{tuAt) 

f— sxn -j- 2 

. . 1 ri 


sinuAt . sin2uAti 

(U ^8 

‘ 4 ■" 32 J 


[— u)At - 2 sin<»JAt + 

sin2tuAt 1 



* w 

••2 

4 •* 




/o^sinu&t cosujAt d(At) = “ /o**^^sin(2iiJAt) d(2uAt) = - ^ [cos (2wAt)]o^^ 

« ~ (1-cos2(»^t) 


/o^sin<^t d(At) ■ /o^^sin<^t d(<*Jdt) “ “ [cos<frAt]o^^ * ^ (1-cosuAt) 


1 

“ /q cos<^t d(oiAt) 


1 , » » “At 

— (sinwAt)o 
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